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ABSTRACT 

^ . We develop an event generator, NLL-QEDPS, based on the QED shower including 

the next-to-leading logarithm correction in the e + e~ annihilation. The shower model 
is the Monte Carlo technique to solve the renormalization group equation so that 
they can calculate contributions of a m \og n {S/m 2 ) for any m and n systematically. 
Here a is the QED coupling, m e is the mass of electron and S is the square of the 
total energy in the e + e~ system. While the previous QEDPS is limited to the leading 
logarithm approximation which includes only contributions of (a\og(S / m 2 e )) n , the 
model developed here contains terms of a(a\og(S / m 2 e )) n , the the next-to-leading 
logarithm correction. The shower model is formulated for the initial radiation in 
the e + e~ annihilation. The generator based on it gives us events with q 2 , which is 
a virtual mass squared of the virtual photon and/or Z-boson, in accuracy of 0.04%, 
except for small q 2 /S. 
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1 Introduction 

In high energy reactions with electron beams, it is important to study radiative 
corrections For this study event generators are indispensable tools. We have 
made the event generator, QEDPS for radiative corrections in the e + e~ anni- 

hilation based on the shower model, which can radiate any number of photons. How- 
ever, this model is limited to the leading logarithm(LL) approximation. In this paper 
we develop a shower model in the next-to-leading logarithm(NLL) approximation 0- 
T0|1 . The magnitude of the NLL order correction is of a 2 /ir 2 \og(S/m 2 e ), which is 



about 0.0001, if y/S is lOOGeV. So the NLL shower might be irrelevant to ac- 
tual measurements. However this is wrong since contributions due to soft photons 
are large. They are estimated to be a 2 / it 2 \og(S / ml) log 2 (E y / y/S) , which is about 
0.005 if the measured energy for a observed photon is lOOMeV. This value is not 
negligible in precise experiments. 

In this paper we limit the event generator in the NLL approximation to the e + e~ 
annihilation, especially to the radiative process on the initial state. Applications of 
the NLL shower to other process such as the Bhabha scattering || are discussed in 
other papers. 

Our study is completely based on the renormalization group equation(RGE), 
which has been developed well in QCD||. First we clarify meanings of the NLL 
order approximation. Let us consider a dimensionless observable F(Q 2 /fi 2 , «o) with 
the mass scale Q 2 and the renormalization point /i 2 . If the coupling ao is small at fi 2 
and the ratio Q 2 /fJ? is large, the RGE shows us that F(Q 2 /fi 2 , «o) can be expanded 
by the coupling constant after summing terms of [cto log(Q 2 / fi 2 )} n for all n. 

F(Q 2 //i 2 ,a ) = F^\a,\og{Q 2 ^ 2 )) + a^ 2 \a,\og{Q 2 /^ 2 )) 

+ a 2 F^(a \og(Q 2 /fi 2 )) + ... (1) 

The first term is the LL order approximation. If the the second term is included, 
the approximation is of the NLL order. 

Sect. 2 contains discussions on the RGE and the formulas needed for later sections. 
In Sect. 3 we apply them to the e + e~ annihilation process and give quantities such 
as the anomalous dimensions explicitly. In Sect. 4 we formulate the shower model 
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briefly. In Sect. 5 we discuss the singular behavior of the NLL order correction. The 
simple perturbative expansion breaks down due to effects of soft photons so that 
more sophisticated techniques are introduced. Then the effective shower model in the 
LL order is given. At this stage the model has the same form as that in the LL order 
model. But some constraints are imposed so that it has included some contributions 
of the NLL correction. Sect. 6 contains results of the effective shower model. In Sect. 7 
we present the explicit form of the second order P function used in the NLL shower 
model. In Sect. 8 we construct the event generator by defining kinematical variables 
in terms of the variable in the shower model. Also some problems on the construction 
are pointed out. The NLL approximation needs the second order coefficient in the 
(3— function of the coupling, but we drops contributions from this coefficient, which 
is discussed in Sect. 9. In Sect. 10 we present a method to compensate results by the 
shower because they contain the Q 2 -independent contribution due to the constraint. 
Sect. 11 is devoted to conclusions and discussions, where numerical results in our 
study are summarized as well as limitations of our model are given. Also we make 
some comments on applications of our model to QCD. 

We present three appendices for some technical parts of our model. Appendix 
A gives us relations between the usual perturbative expansion and the RGE. In 
Appendix B we discuss the approximation that is made in order to get analytical 
expressions for results by the shower model. In Appendix C we present compact 
descriptions for the shower algorithm, which we apply to compensating results for 
the Q 2 -independent contribution. 



2 Renormalization group equation 

In the RGE0, the value of the coupling depends on /x 2 , so that it is not constant, 
but the function of /i 2 , a(/i 2 ). Then the derivative of the coupling by fi 2 is given by 
the function of the coupling only. 




3 



Solving this equation, we obtain that 

log(Q 2 //i 2 ) = f J^da, (2) 

a = a(Q 2 ), a = a(fx 2 ). 

Here a is called the running coupling. If one applies the RGE to the dimensionless 
observable F(Q 2 //i 2 , a ), the following equation is obtained. 

(fi 2 d/ dfi 2 + P(a )d/da - 7 (a ))i 71 (Q 2 /V, «o) = 0. (3) 

7(a) is the anomalous dimension which is the function of the coupling a only and 
depends on a process. Then we solve this equation to obtain that 

F{Q 2 /fi 2 ,a ) = exp(- [" i{a)/P{a)da)F{l,a). (4) 

In the NLL order, 

(3(a) = p t a 2 + (3 2 a 3 , 
Pi = 77-, # ' 



3vr' ^ 2?r 2: 
7(a) = 7i« + 7 2 a 2 . 



Solving Eq.(@) on the running coupling, 

A log(g 2 /^ 2 ) = "(---)- M«/«o). (5) 
a cto pi 

If a is small but a log(Q 2 //i 2 ) is large, we keep any term of [<a \og{Q 2 / ji 2 )} n and 
drop terms of a§ \og(Q 2 / fi 2 )(K > 3). Then we obtain the explicit formula for the 
coupling at Q 2 in the NLL order. 

^0 a (3 2 logjl - atoPi \og(Q 2 / jJ 2 )} ( ) 

a 1 - aoA log(g 2 //" 2 ) /Si 1 - aoA log(Q 2 //i 2 ) 1 J 

The integral inside the exponential of (f|) on the anomalous dimension is carried 

out in the NLL order 



7i 1 /- / \ , ,72 7i# 



J 7 = / j(a)//3(a)da = -^log(a/a ) + (7^ - -^)(a-a ). (7) 

Ja Pi Pi Pi 

The Q 2 dependence of F(Q 2 / /1 2 , a ) can be expressed by p\, p2, 71, 72 and notations 
used in Appendix A. 

F(Q 2 //, 2 ,ao) = exp[- f" da^]F(l,a) 

Jan Pi a) 



r 7l , /«x ,72 7l#2w_ v,/, . ,Q _x 

exp[-— log(— ) - (- _)( a _ ao ) (/ + / a ). 

Pi a Pi pi 



To obtain the formula for the Q 2 dependence, we take the ratio of F(Q 2 /fi 2 , ao) and 
F(l,a ). 

F(Q 2 //j 2 ,a )/F(l,a ) 

= expi-f log(^) -(f-l±fo)(a- aoWo + #*)/(/„ + fla*). (9) 
Pi "o Pi Pi 

Here since (fo + f®a)/(fo + /fao) can be approximated by 1 + ^(a — ao), 

2 I ..2 



F(g7 / i i ,a )/F(l,a ) 

exp[-— log(— ) - (- ^-)(a - a )][l + —(a - a )]. (10) 

Pi «o Pi Pi Jo 



Then we expand terms of a— a in the exponent and drop terms of (a—a ) K (K > 2). 

F(Q7/i 2 ,ao)/F(l,a ) 
= exp[-^log(-)]{l + (-^ + ^ + ^)(a-a )}. (11) 



Pi a Pi Pi /( 
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Finally we describe the explicit expression of the Q -dependence for F(Q //i , ao). 

F(Q 2 //i 2 ,a )/F(l,a ) 

r 7i i fll (n 2, 2\\ . 7i/?2a log(l -a /?ilog(Q 2 //i 2 )) 

= exp[— log(l - a pi log(Q /// )) + 



L Pi taV ^ ^ i_ ao ^bg(gV/x 2 ) 

72 7i& A° PiallogjQ 2 / n 2 ) 
Pi PI fJl-aoPrlogiQ 2 /^ 2 ) 



The above is the fundamental equation for the NLL shower model. 



3 Annihilation 

In this section we present explicit formulas for the annihilation cross section. 
When o"o(Q 2 ) is the bare cross section, i.e. the cross section without any radiative 
correction, the observed cross section is expressed by the structure function D(x, Q 2 ) 
and the coefficient function C(x,a) JTT| |. 



da obs (S,Q 2 ) 



a (Q 2 )^ ^ d Xl — f 1 dx 2 —D(x 1 ,Q 2 )D(x 2 ,Q 2 )C(z,a), (13) 
aQ z b Jo X\ Jo X2 



T Q 2 Q 2 q 

-, T = — , S = X1X2O, 



X\%2 s S 
where S is the total energy squared. 

In oder to solve the RGE in the NLL order analytically, we take moments of 

Eq.©. 

da obs (n,Q 2 ) cr (Q 2 ^ 



dQ 2 Q 2 
Here they are given by taking moments 



D(n,Q 2 )D(n,Q 2 )C(n,a). (14) 



da obs (n) r 1 da obs (S,Q 2 ) 

~ ' CLT — — r , (15; 



dQ 2 Jo dQ 2 

D(n,Q 2 )= / —x n D(x, Q 2 ), (16) 
Jo x 

fi dz 

C(n,a)= / — z n C(z,a). (17) 
Jo z 

D(n,Q 2 ) corresponds to the exponential term in Eq.(|J), while C(l,a) does to 
F(l, a) so that 



D(n,Q ) = exp - / da], 

Jan P(a) 



>ao (3(a) 

7(71, a) = 71(71) a + 7 2 (n)a 2 , (19) 



C(n,a) = l + ^C 1 {n). (20) 

Z7T 



Then we obtain the Q 2 -dependence of the cross section, as in Eq.(|T 

:<2 — 37^ — = ex Ph 2 — ^iog(— )] 



x {l + (-2^ + 2i^ + ^ 1 H)(a-a )}. (21) 



<7o(Q 2 r dQ 2 fa & W 

, 72(71) 7i(n)/3 2 

A A 2 ' 2tt' 

Here we summarize 71 (n), 72 (n) and other quantities needed for the annihilation 

process. In order to get the DGLAP equation, which is used in the shower model, 

we replace the variable fi 2 by Q 2 in the RGE, so that the equation of D(n,Q 2 ) on 

Q 2 is 

Q 2 ^D(n, Q 2 ) = (Q 2 ^ + p(a)-^)D(n, Q 2 ) = - 7 (n, a)D(n, Q 2 ). (22) 



Of course D(n,Q 2 ) of Eq.(|T8|) satisfies the above equation. The inverse Mellin 
transformation of the above is called the DGLAP equation. 

Q 2 dD( f n f ] = f ^P(x/y, a))D(y, Q 2 ). (23) 

Here P(x, a) is called P function, which is defined by the anomalous dimension 

7(n,a); 

j(n, a) = — dxx n ~ 1 P(x,a). 
Jo 

P(x, a) is divided into the LL and the NLL terms. 

P(x, a) = ^P^(x) + (^) 2 P (2) (x), (24) 
= P + ( x ), 
P(x) = • 

1 — X 

Here we introduce the +notation for a function f(x). 

f+{x) = f(x)-6(l-x) [ 1 f(y)dy. 

Jo 

The explicit form of 71 (n) is useful for examining the model. 

7l (n) = te-'P«(,) = l(2S l( „ -D-1 + 1 + -L.). 

Here we introduce functions on the summation. 

n y 

Sm(n) = ]T — for n - X > 

k=l K 

S m {n) = for n < 0. 

The NLL term in the moment expression is not so compact so that we present only 
P( 2 )(x), which is found in Refs.[0],|T§g. 

l2(n) = -\- 2 f'dxx^P^ix), 

P {2 \x) = P a+ (x) + P b+ (x), 

3 

P(x) a = — P(x) log(a;) log(l — x) — ( h 2x) log(x) 

1 — x 

- i(l + x)log 2 (x)-5(l-x), (25) 



1 We drop terms of branching of e into e + , which are quite small. 



P b (x) = ~[P(x)(- log(x) - 5) - 2(1 - x)}, (26) 

C(a;,a) = <J(l-x) + ^Ci(x). 

Ztc 

Here C*i(x) has been calculated for the Drell-Yan process in QCD[jll], which is the 
sum of C ep+ in the deep-inelastic scattering and an additional term AC(x). 

C 1 (x) = 2C ep+ (x)+AC(x), 
C ep (x) = P(x) [log(^) - | + (27) 

A(z) = 2P(x) log(l - x) - — — - 6 - 4x, (28) 

AC(x) = A + (x) + 6(1 - x){-- + ^). (29) 

In our model we make use of the scheme dependence, which says that only the 
combination of — 2 :!2 ^ + ^Ci(n) can be predicted by the RGE, but each quantity 
is not. By making use of this freedom, we can put C\(n) = 0. 

-^ + £*W-^ + £ft W -^. (30) 



That is 



Or 



7a(n) = 7i(n) - ^d(n). (31) 



P( 2 ) ( X ) = p(2) ( X ) + yr^d (x) . (32) 



In our model we use P^(x), which means that the hard cross section is not em- 
ployed. 



4 Shower model 
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The shower model that we call in this paper is the Monte Carlo method to solve 
the DGLAP equation by repeating the branching that the electron branches out into 
the child electron and the photon, where a variable x and a virtual mass squared K 2 
are generated. Here x is the energy fraction of the child to one of the parent, while 
K 2 is the absolute value of the virtual mass squared of the child. So the moment 
distribution of x^ 1 = (x±X2 ■ ■ • xl)™ -1 calculated by the shower model agrees with 
the analytic result of the RGE completely within the statistical error. Here 
at the i— th branching and L is a number of branchings in one branching process. 

The shower model needs the scheme to cutoff the infrared singularity, though it 
is arbitrary. Since we would like to apply the shower model to the event generator 
that produces electrons, photons and other particles in simulations, we adopt the 
following cutoff scheme. 

x < 1 - n 2 /K 2 . 

The definition of x, which is necessary to construct the generator, is given in Sect. 8. 

Also our shower model employs the double cascade scheme, in order that the 
electron and the positron make the branching process independently In 
this scheme we impose the constraint to x. 

l-x> K 2 /Q 2 . 

By using these constraints, we can apply the Monte Carlo method to generate x 
and K 2 , as described in Ref.p|. 



5 The singular behavior of the NLL order correction 

In Sect. 3 we presented the P function in Eq . (|32|) , which we use in the shower 
model. As x — > 1, the most singular behavior of this function is 
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This is dangerous, since it can become larger than the singular LL term, which is 

a(K 2 ) 2 
2tt l-x' 

In order to make the branching stable, it was suggested to use the running coupling 



a 



((1 —x)K 2 ) instead of a(K 2 ) in order to include term of a n ° g 1 _ x into the model 
in Ref. [0]. By this replacement, the singular behavior of P function is 

a(K 2 ) 2 + a(K 2 )^log(l-x) 



2tt 1 - i 2tt 1 — x 

«o 2 t a(K 2 )^ 2 \og(l -x) 



2tt[1 - PxOq \og(K 2 / /j, 2 )] l-x v 2tt ' l-x 

2 +{A _ 47r(3i){ ^K 2 ), 2 lo g (l-x) 



2n[l- (3 1 a \og((l-x)K 2 /i2 2 )]l-x v ^ /v 2tt ' l-x 

+0{al). 

In the deep inelastic scattering, indeed A = 4n/3i, so that the dangerous term 
disappears after the replacement. On the other hand we have A = 87r/3i in the an- 
nihilation. By this method 47r/3i of them can be included into the running coupling. 
We have to include the remnant 4n/3i terms into the effective LL form in order to 
remove the dangerous term in the second order P function. This can be done by 
taking account of the kinematical constraint in the annihilation as follows. 
In the annihilation, both electron and positron to radiate photons. 

e-(Px) + e+(P 2 ) -> e-( Pl ) + e + (p 2 ) + X - 7 (g) + X. 

This implies that the spacelike virtual electron(pi) and positron (p 2 ) annihilate into 
the virtual photon(g). Although we can calculate q 2 by pi, p 2 , as seen in Fig.l, we 
make the approximation that 

q 2 «z w x 62 (l-ti)(l-t 2 )S. (33) 

Here t\ = —p\IQ 2 {ti = —pl/Q 2 )- The accuracy of this approximation is discussed 
in Sect. 8. 

Therefore our shower model calculates moments with respect to the variable 
Xb(l — t) . The n— th moment (xb(l—t)) n ^ 1 should agree with the structure function 
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D s {n,Q 2 ) on the electron(See ref.Q]' 

dt 



D s (n,Q 2 )=U(Q 2 ,fi 2 ) + [ -(l-t) n - l n(Q 2 ,tQ 2 ) 

Jo t 

x f 1 dx a ^ 1 ~ %)tQ2) P(x)x n - 1 6(l -x- t)0((l - x)t - e) 
Jo 2n 

r /•*£# f, ,Qt((l - x')t'Q 2 ) , , ln _ x 

x exp{ / — / da; — P{x){x - 1) 



/o £' J 2tt 
x^(l-x'-t')^((l-^)^-e)}, (34) 



n(g 2 ,gg) = exp [- / -/da; 

J olio 2 t J 



dt f a((l — x)tQ 



2a 



xP(i)fl(l-a-t)0((l-x)t-e)]. (35) 

Here e = fi 2 /Q 2 . U(Q 2 ,Ql) is the non-branching probability that the electron does 
not branch for possible virtual mass squared between Q 2 and Qq. In Eq. (|34]) the first 
term represents the no-branching case so that the moment is unity for any n. The 
front term on the exponential does the last branching, while the exponential appears 
after repeating branchings. The reason for the special form on the last branching is 
that there the virtual mass squared is involved in the moment, as seen in Eq . (j33|) . 

In order to obtain the expression that is possible to be calculated analytically, 
we approximate (1 — t) n by 6{l/n — t). 

- 1 dt f 1 a((l-x)tQ 2 ) 



D s (n,Q 2 ) = exp{ / — [ dx 
Jo t Jo 



2tt 

re-l 



xP(x)[(x(l - t)) 71 - 1 - 1]0(1 - x - t)0((l - x)t - e)}. (36) 

The error due to this approximation is discussed in Appendix B. We write D s (n, Q 2 ) = 
exp[I s (n, Q 2 )], and we perform the integrals for I s (n, Q 2 ). 

L(n,Q 2 ) = 

= f 1 dz{z n - x - 1) f 1 - f 1 dx5(z - x(l - t) f^ l ~ x)tQ2) P{x) 
Jo Jo t Jo 2tt 

x [(xil-t)) 11 - 1 -l]6(l-x-t)6((l-x)t-e). (37) 

Then we integrate I s (n,Q 2 ) over x and t, where we assume that n is not so large 
so that contributions from regions of z ~ 1 can be neglected. Also we neglect 0(e) 
terms. Finally we obtain 

h(n,Q 2 ) 
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= -L f 1 dx(x n - 1 - l)P(x)±- log(— ) + ^ rfx^ 1 - l)AP(x), (38) 
An Jo pi ao Ztc Jo 

where 

AP(x) = 2P(x) log(l - x) - p ^ log a; 

+ (l + x)log(l + Vx) + -(l + x)logx- v^+a;- (39) 

In Eq.([38|), the first term is the LL order while the the second order is the NLL order, 
though it contains the Q 2 -independent contribution. By noting that a = (a— ao)+ao 
and 

f %(tQ 2 )] 2 = Ua-a ), (40) 

one can see that the second term has (a/2n) 2 log(l — x)/(l — x) and its coefficient 
is 87r/?i as expected. 



6 Results by the shower in the effective LL order 

Summarizing the discussion in the previous section, Eq.fl3"BD was derived by adopt- 
ing the three schemes: 

1) to use (1 — x)K 2 for the argument of the running coupling, i.e. at the branching 
we employ the following coupling, 

a((l - x)K 2 ) = — 

U ) ) l-/3ia log((l-x)^ 2 //i 2 )' 

As pointed out in Ref. [0, (1 — x)K 2 is about the transverse momentum squared at 
the branching. 

2) the double cascade scheme ||], where we impose the constraint of 

(l-x)>K 2 /Q 2 . 



3) to define q 2 in the annihilation as 

q 2 = x bl (l - h)x b 2(l - t 2 )S, 
12 



Xj, — 1 X X\ X • • • X Xl- 



The moment of Xb(l — t) by this shower in Eg. (p8|) is calculated in the analytic 
form. 

Is(n, Q 2 ) = J- log(^) [-2^(n - 1) - i - -L- + |] 

+ ^[2(S 2 (n - 1) + S 2 (n - 1) + iSxCn) + -^^(n + 1) - ~) 
27r n n + 1 4 

- S 2 (n - 1) + -(-Si(n) + 5i(2n)) + ——(-S^n + 1) + S'i(2n + 2)) 
n n+1 

1,J_ 1 2 1 

2^n 2 + (n + l)^ 2n + l + n + l J ' 

(41) 

Some comparisons between analytic calculations, Eq.([|l|), and results by the shower 
model are shown in Table 1. There we assumed that fj 2 = 0.25 x 10 _6 GeV 2 and cto = 
1/137. In Monte Carlo simulations a total number of events is 10 8 and errors are 
estimated by calculating the variance in 10 data sets of 10 7 events. The agreement 
of order of 10 -5 justifies our discussion. 

The structure function given by Eq.fl4TD has the Q 2 -independent contribution, 
Df(ri), which does not vanish at Q 2 = jj 2 because the term of a in Eq.(|41|) remains 
then. Df(n) is obtained by replacing a by ao in Eq. (|38|) . 

DHn) = exp{— f 1 dx{x n ~ x - l)AP(x)}. (42) 

27T JO 

In order to calculate the absolute value as well as the Q 2 -dependence for the structure 
function by the shower model, we have to compensate it for the contribution due to 
Df(n). The method for the compensation is discussed in Sect. 10. 



7 Effective P^(x) 

In this section we give the second order P function p( 2 > e fl \x), which is used in 
the shower model. Since our model imposes that the coefficient function is zero 
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so that the NLL contribution can be given by P^ 2 \x) in Eq.(^). As discussed in 
Sect. 5 and 6, our shower model contains the NLL contribution, AP(x), through 
a((l — x)K 2 ), the double cascade scheme and the definition of q 2 , as described in 
Eq. (PS]) . In the NLL shower we employ these schemes so that the structure function 
D(n, Q 2 ) = exp[I(n, Q 2 )} is given by 

r 1 r 1 dt 

J(n,Q 2 )=/ dx -^(l-x-t)^((l-x)t-e)[(l-t) n - 1 x n - 1 -l] 
Jo Jo t 

{ ^{l-x)Q 2 ) p{x) + a(t(l - x)Q 2 ) 2 p{2) eff {x) y (43 ) 
2n (2n) 2 

By performing the integral similar to one in the previous section, we have 

I(n, Q 2 ) = log(-) f dx{x n ~ l - l)P{x) 
2-k/ji a Jo 



+f /V"" 1 - l)AP(z) + t^^(« - « ) f dx{x n ~ l - l)P^ 
2n Jo (2n) z pi Jo 

-log(— ) C dxix 71 " 1 - l)P(x) 
1 «o Jo 

f 1 dx^x 11 ' 1 - l)[P ( V eff (x) + 2n(3 l AP(x)} 
Jo 



eff( x ) 



2np 

fl 



+ ^ / dx(x n ~ L - l)AP(x). (44) 

Z7T JO 

Here we neglected terms of order of or a 2 . P^^x) in Eq. (|32D should equal to 
P^ eff (x) + 2n(3 1 AP(x) . Therefore P^ eff (x) is given by 

p(*hff( x ) = p(2)^ _ 27r p 1 AP(x) 

= P a (x) + P b (x) 

+ 2tt/3i[-— — - P(x) log(x) + - (1 + x) log(x + y/x) + - x],(45) 

1 — x 1 — X 

which is free from the singular term log(l — x)/(l — x), as is expected. 



8 Event generator 

In this section we present the event generator based on the shower model, which 
was described in Sect. 4-7. In order to determine four momenta of the produced 
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particles we define i to be a fraction of +(— ) component of lightcone variables for 
electrons(positron) . 

x = pr(p-)> (46) 

P± = -fi- (47) 

Here four momenta are denoted as (p x ,p y ,p z , E) and P_(P + ) denotes a lightcone 
variable of the initial electron(positron) f2|. At the branching of e~(y, — K 2 ) — > 
e~(xy, —K' 2 ) + 7(y(l — x),0), the momentum conservation imposes the following 
equation. 

_ * = ~JS1 + _f_ (48) 
x x(l — X) 

where kr is the transverse momentum, {p x ,Py)- Here note that the electron during 
the branching process is spacelike. Our cutoff scheme, x < 1 — fi 2 /K' 2 , equals to 
k T 2 >^ 2 , if K 2 ^K' 2 . 

Using an arbitrary azimuthal angle 0, K 2 ,K' 2 and x one can determine four 
momenta of the electron and the photon after the branching. 

_ K' 2 4- k 2 -, 

p e ,, = (xP + , ±^,hr), (49) 

x 

Vlil = ((l-x)P + , r ^,-A r T ), (50) 

k x = k T cos(p, k y = k T sin (ft, (51) 

k 2 T = (l-x)(-xK 2 + K' 2 ). (52) 

These equations determine four momenta of all particles completely. This implies 
that in the annihilation process the four momentum q of the virtual photon and/or 
Z-boson is the sum of momenta p\ and P2 of the electron and positron after the 
branching process, that is q = p\ + p 2 (See Fig.l). Then the virtual mass squared of 
the four momentum is 

q 2 = pl+pj^ 2 Pl p 2 

= -K 2 - K 2 + 2(p 1+ p 2 - + P1-P2+ - PitP2t) 

= -K\ - K\ + 2(x bl P l+ x b2 P 2 _ + ~f + / 1T ~f + /* T - 2p lT p 2T ) 

2x bl P 1+ 2x b2 P 2 - 
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The variable r = q 2 /S equals 



T 



t\ — t 2 + X b \X b2 + 



x bl x b2 S 2 



(54) 



Here we used t = K 2 /S. In the generator the ratio r is given by the above equation, 
but not by x b \x b2 . Although the RGE predicts moments on r, as described in Sect. 3, 
r of ([54]) is not a good variable for the shower model, because it gives us moments on 
T ' = x bl {l — ti)x b2 (l —t 2 ), as discussed in Sect. 5. We present the detailed discussion 
on the accuracy of the generator, which uses r. 

Let us discuss differences between moments on r and r'. First note that the last 
term, 2pi T p 2T , in the above equation is zero if averages are took, because angles 
between these vectors are arbitrary. Next r can be negative while r' is always 
positive. Of course the negative r is unphysical so that the event with the negative 
t is abandoned. We introduce a variable r that can be negative and whose moments 
are possible to be calculated analytically, which is 



In this definition r might be negative. A case of being negative is counted as an 
event, but these negative values are replaced by zero. Also the generator can fail to 
make four momenta for the virtual photon and/or Z-boson, i.e. r is negative. The 
failed case is counted as an event and r is set to be zero. Results on moments of 
r, r and r' for a total number of 10 8 are presented in Table 2. Differences on data 
for t in the generator and r in the shower are quite small and less than 10~ 5 . So 
we can conclude that moments of r in the generator is accounted by those of r. We 
can estimate analytically differences of moments on r and r', which are of order of 
10~ 3 and decrease rapidly as n rises. They are 



T = (x b i - ti)(x b2 ~ t 2 ). 



(55) 



1 r 1 dt r 1 
D d (n,Q 2 )=exp{2— - dx6{t(l-x) 

Z7T JO t JO 

a(t(l - x)Q 2 )P(x)[(l - t)™" 1 ^- 1 - 1]} 



6)0(1 



X — t) 




a(t(l - x)Q 2 )P(x)[6{x - t)(x - t) 



!]}• 



(56) 



1G 



Here for simplicity we neglect P^ e ^(x). If we neglect the running effect on a and 
approximate Dd{n, Q 2 ) by the difference between the first terms in expansions by a, 

D d (n, Q 2 ) « ^ [' % t dxP(x)d(t(l -x)- e)6((l - x) - t) 

Z7T JO t JO 

x [^(l - ty- 1 - 9(x -t)(x- t)"- 1 ]. (57) 

Under these approximations, Dd(n,Q 2 ) is independent of Q 2 . The values of the 
expression Q57|) presented in Table 3 can account for the difference on the moment 
of r' and r. For an example, the moment of n = 2 in Table 2 is 0.93212, and 0.93118, 
while Table 3 shows the difference 0.00093, which agrees with the difference. 



9 /3-function 

In this section we estimate contributions by the second order correction /3 2 m the 
beta function. Using Eq.(|) for the running coupling, a — ao is given by 

ftaglogWV) 

° " a ° = nmS(W)' (58) 

Inserting the explicit expression of the running coupling, the integral on the anoma- 
lous dimension 7 7 in Eq.(|7]) becomes 

T 71 i 1 



A b i-«oAlog(Q 2 //i 2 ) 



a /? 2 7i log[l - apfh \og(Q 2 /n 2 )} 7 2 jifo gia|log(gV) 
/?i l-a ftlog(Q 2 M 2 ) Pi A 2 J 1 - a 0/ 9i log(Q 2 /> 2 ) ' J 

In the QED process ao/?i log(Q 2 / /i 2 ) is small so that the logarithm of the second 

term is approximated as follows. 

log[l - otoft log(Q 2 //i 2 )] « -«oA log(Q 2 //^ 2 )- (60) 

This leads to cancellation of the terms with /? 2 . 

T ,7i, 1 2 MQV) rfi1 x 

7 ~ ft ° g 1 - a /?i log(Q 2 //i 2 ) + 72C " 1 - aoft log(Q 2 //i 2 ) ' 1 j 
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The leading term in the neglected terms is 

lifa^l MQ7/X 2 ) 
2 l-a /3ilog(Q 2 //i 2 )' 

which is less than 10~ 6 if one uses actual values for oto,fa,fa, 7i and Q 2 /fJ> 2 = 
(100GeV/0.5MeV) 2 = 2.5 x 10 11 . Therefore we can neglect the term with fa safely. 



10 Q 2 -independent contribution 

As discussed in Sect. 5, the shower model contains the Q 2 -independent contribu- 
tion. In other words, the structure function is not 8(1 — x) at Q 2 = fi 2 , but its 
moment is given by D'(n) of Eq.(f42|). In study of the radiative corrections on QED, 
the absolute value of the cross section, not its Q 2 -dependence, has to be calculated. 
So we would like to compensate the structure function for the the Q 2 -independent 
contribution D^(n). On the moment we only divide results by D'(ri). 

D cmp (n,Q 2 ) = D(n,Q 2 )/D f (n). 

Here D cmp (x, Q 2 ) is the structure function after the compensation. However, in the 
generator we need the inverse Mellin transformation. The product of moments is 
equivalent to the convolution integral of the function in the transformation so 

D cmp (x,Q 2 )= f 1 ^D(y,Q 2 )D^(x/y), 
Jx y 

(D f (n)Y l = [ dxx n - 1 D 7 (x). (62) 
Jo 

In the shower the convolution integral is realized by the procedure that x of the 
initial electron is fixed according to the probability D*(x) and then we make the 
branching process, which induces the structure function. 

For performing this we are confronted with three problems. First the explicit 
form of D*(x) is difficult to calculate. Here we use the shower algorithm to get x, 
which is described in detail in Appendix C. 
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Second the function D f (x) is not positive necessarily so that it could be used as 
the probability. Equivalently the splitting function —AP(x) in 

(D f (n))- 1 = exp[-^ j 1 dx{x n ~ l - l)AP(x)], 
2tt Jo 

is negative at x near zero. Since —AP(x) is concentrated near x ~ 1, we must 
be contented with an approximated AP A (x), which is modified near x ~ 0. The 
approximated AP A (x) is given by 

AP A (x) = 6{x - x c )AP{x). (63) 

Here x c is fixed to satisfy the condition, 

dyAP(y) = 0. 



f 
Jo 



The third problem is that the total energy squared S of e + e system is changed 

f f f f 

by x{x 2 S, where x{ and x 2 are x— fractions of the initial electron and positron 

after the compensation by D^{x). Since errors due to this problem are very small 

numerically, we neglect this problem. 

Next we examine numerical results related to D^(n). First we compare (D* (n)) -1 

with 

D f ' A (n) = exp[-^ [\x n ~ l - l)AP A (x)}, (64) 
2it Jo 

in Table 4. One finds the agreement between results by the shower algorithm and the 
moment (D^n))' 1 except those for small n. The differences are less than 0.01%. In 
order to confirm that these differences are due to the second problem, we calculate 
the moment D^ A {n) through the numerical integration to show these results in the 
same table. The agreement between values of the second and the third column 
supports our discussion strongly. 



11 Conclusions and discussions 
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In this paper we have formulated the shower model including the NLL correction 
in the e + e~ annihilation and developed the generator. Results on the Q 2 - dependence 
of moments by the generator are summarized in Tables 5, 6 and 7. Here the moment 
for n — 1 is fixed to be unity because it is the event normalization, so that the 
analytic value is the ratio between the n— th moment and the first moment. Analytic 
results in the NLL order are 0.03% of the LL order ones for small n, but they 
increase as n rises and are about 0.2% for n = 100. The effect of the NLL order 
is small, but could not be neglected in precise experiments as at LEP or in future 
colliders. Results of our shower model agree with the analytic calculations of the 
NLL oder in accuracy of 0.04%. The agreement in the generator is worse for small 
n, because events might not satisfy with the kinematics constructed by the shower 
model. Simply speaking, the constructed value for q 2 becomes negative in these 
events. Table 7 shows differences between moments by the shower model and the 
analytic one, and those between moments by the generator and the analytic one in 
the NLL order. Values in the table indicate the magnitude of the systematic error 
in the present shower and the generator in the NLL order. 

Next we mention some comments on limitations of our model. First the accuracy 
found in q 2 distributions may not be common to other distributions such as the 
transverse momentum distributions of radiated photons. One of reasons is that we 
do not include the cross section for the emission of a photon with a large transverse 
momentum. 

Also we neglect effects by the three-body decay in the shower because there is 
a less interest on detailed distributions of photons. Further we neglect the mixing 
p( 2 \x), which is a contribution that the electron radiates into the spacelike positron 
with the pair creation, since its effect is expected to be very small. 

Third our shower model is limited to the non-singlet case where there are no 
contributions by radiations with the spacelike photon. A reason for this limitation 
is that in experiments one can exclude events with the electron positron pair easily, 
which correspond to pure singlet radiations, as well as that they are quite small. 
Finally notice that we use S for the mass scale of the RGE, but not q 2 . 

In this study we have examined our model in the moment form. But from a 
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experimental view, analyses in the x-space are desired, which will be discussed in 
coming papers. Also we have to discuss the accuracy of our generator in detail by 
applying it to several realistic processes such as the muon pair production or the 
Z— Higgs production. 

Finally we would like to stress that our study on the NLL shower is quite impor- 
tant for QCD, where the NLL shower has been developed. Because there has been 
no study on the (^-independent contributions by the shower algorithm in QCD. 
Also precise discussions have led us to the deeper understanding of the shower mod- 
els. Therefore our study stimulates the interest on further developments of QCD 
showers. 
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Appendix A 

In this appendix we explain results of the RGE in terms of the perturbative 
expansion. In the perturbative expansion by the coupling a at /x 2 , the dimensionless 
quantity F(Q 2 //i 2 , a ) is calculated up to some order of a . 

F(Q 2 //i 2 , oo) =f + fi(Q 2 /n 2 )a + / 2 (Q 2 //i 2 )«o + • ■ • • (65) 

Here we assume that the mass scale is only Q 2 . In order to compare results of the 
RGE with those of the perturbative calculation, the ratio F(Q 2 /[i 2 , a Q )/F(l, ato) 
should be used, because the RGE can calculate the Q 2 - dependence only. The ratio 
is given by 

F(Q 2 /fi 2 , a )/F(l, a ) = l + [(fi(Q 2 /» 2 ) - /i(l))//o]« + 

UiQ 2 /^) - / 2 (l))//o + /i(l) 2 //o 2 " /i(tfV)/i(l)//?]«8 + " " " • (66) 

The physical quantity contains only log(Q 2 //i 2 ) and the order of power of the loga- 
rithm is less than the order of the coupling constant. Therefore 

fi(Q 2 /^) = fi+f^og(Q 2 /^ 2 ), (67) 
/ 2 (Q 2 //i 2 ) = / 2 ° + fl MQV/i 2 ) + fl log 2 (g 2 //i 2 )- (68) 
If the ratio F(Q 2 //i 2 , a )/F(l, a ) is expressed in terms of /°, fl, f®, fl and /f , 

F(Q 2 //i 2 , a )/F(l, Qt ) = 1 + aj-j- log(Q 2 //i 2 ) 

Jo 

+ log(Q 2 //, 2 ) + & log 2 (Q 2 //, 2 ) - log(g 2 //, 2 )}. (69) 

Jo Jo UoJ 

The result by the RGE has been discussed in Sect. 3 and is summarized by Eq. (p^) . 
If we drop terms of (n > 3) in this equation, 

F(Q 2 /fi 2 ,a )/F(l,a ) 
« 1 - a 7ilog(g 2 //i 2 ) + a 2 log 2 (g 2 //i 2 )[i( 7 i) 2 - 7 i/3i] 

+ a 2 log(Q 2 //i 2 )[- 72 + 4^]- (70) 

Jo 



Since the perturbative expansion (|69"D should agree with (|70[) by the RGE, we obtain 
relations between 7, f3 and /. 

fi 

-j- = -7i, 
Jo 
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f = ^(7i) 2 -7iA, 
Jo 2 

li_f±IL _ A°A 

/o (/o) 2 " 7 fo ' 
If we neglect the running effect i.e. /3i = 0, we have the simple relation for f\ 

and /|. 

/ 2 2 = (71) 



Appendix B 

In this appendix we discuss the analytic expression for the moment of Xf,(l — t) 
in the shower model. We need some approximations in order to get the analytic 
expression for the moment, which is defined in Eq.(^) in Sect. 5. The approximated, 
but analytic expression is given by the expression from which we obtain Eq.(^). 
A conclusion in this appendix is that the difference between the expression and 
(j36l) is of order of a 2 and of n°, and does not have the Q 2 -dependence except one 
due to the running effect. 

First note that the expression (|36D equals to 

n(QV) + n(QV) f'ddeMt^ t dx a{{l ~ x)t ' Q2) P(x) 

Jo at Jo t Jo Zti 
(x(l - l/)) n ~ l 0(l -x- t')9((l - x)t' - e)} 

Jo t Jo Zn 

(x(l - t)) n - l 6(l -x- t)9((l - x)t - e) exp{ f — /* dx' ^ ~ X '^' Q ^ P ( x >) 

Jo v Jo 2n 

{x'{l - 0)^(1 -x' ~ t')0{{l ~ x')t' - e)}. (72) 
Since the factorization can apply to the non-branching probability, we have 

Tl(Q 2 ^ 2 )=Tl(Q 2 ,tQ 2 )Tl(tQ 2 ,fi 2 ). 
Then the expression (^) becomes 

n(QV)+ f yn(Q 2 ,tg 2 ) f dx a{[l ~ x)m p(x) 

Jo t Jo Z7T 
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(x(l - t)) n - l 6{\ -x- t)9((l - x)t - e) exp{ f — C dx' ^ P(x') 

Jo v Jo 2ix 

[(x'(l - t')) n - 1 - 1}9(1 - x 1 - t')9({l - x')t' - e)}. (73) 
We examine a difference Dif f = (|3^) -(|36|). 
Diff 

= f 1 -n(Q 2 , tQ 2 ) f 1 dx a ^ l ~ x)tQ2) P(x)(x(l - -x- t)9((l - x)t - e) 

Jo t Jo 2ir 

{ exp [ f % f 1 dx' a ^ l ~ X ' )tlQ2) P{x'){x' n - 1 - 1)9(1 - x' - t')9((l - x')t' - e)] 
jo t Jo 2n 

r ft dt' f 1 , ,a((l - x')t'Q 2 ) „. .... i 
- exp/ — / ds' U ; V ^ ^ l-Or 1 -! 

JO £ JO Z7T 

^(l-z'-OW-^K-e)]}- ( 74 ) 
If the exponential is expanded by a, 

Diff 

~ r 1 ^ r 1 ^ ^(l-a;)^ 2 ) _ t))n -i B ^ _ x _ t)e ^ x _ x)t _ e) 

Jo £ Jo 2n 

[f^j 1 dx M1 ~ x ' )fQ2) p(x')x-\i - (i - tr- 1 )) 



o f Jo 271 

9(1 -x'-t')9{{\ -x')t'-e)}. (75) 

Here we would like to show that this difference is not proportional to log(e) and 
finite as n increases. But it is difficult to obtain the analytic expression for Eq.([75D 
so that we calculate the simpler expression by replacing the running coupling by 
otQ. This approximation could not change the essential property of Eq . (|75|) . If this 
expression is finite as e becomes zero, it is not proportional to log(e). So we set e 
zero in Eq . (]75f) . Then the difference is approximated by 



Diff*(^) 2 t t dxP(x)(x(l-t)T- 1 
Ztt Jo t Jo 

d df r 1 -* 



x / ™ / dx'P(x')x' n - 1 (l - (1 - t') n - 1 ). 
Jo t Jo 

(76) 



This integral is possible to be expressed analytically, which implies that Eq.flTop is 
finite as e goes to zero. But this expression is too long to understand the property. 
In order to confirm that Eq. ( [To] ) is finite as n increases, we present numerical values 
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of them for various n in Table 8, which supports our conclusions. Summarizing this 
appendix, the error by the approximation on D s (n,Q 2 ) is of order of a 2 and n° so 
that we can neglect this difference safely. 



Appendix C 

In this appendix we present the shower algorithm to generate x according to that 
the moment is given by Eq. flS^) , which is 

exp[ f 1 dx{x n ~ l - l)I f (x)l I f (x) = -^AP A (x). 
Jo ' 2ir 

We introduce D (n, x max ), which is 

D J {n,x max ) =exp[ dxx^F {x)). (77) 
Jo 

Since the exponential is expanded to be the infinite series, 

D f (n, x max ) = exp( f*™ dxx^I^x)) = £ [ f dxx n ~ l P {x)] k / k\ 
Jo k=0 Jo 

°EL fXmax rxi rx k -i 

= V / dx 1 xr 1 I f (x 1 ) / dx 2 x n 2 - 1 I f (x 2 ) ■ ■ / dx k x n k - l I(x k ), (78) 
£r Jo Jo Jo 

we can separate contributions of the no-branching and the branching with one or 
more particles. 

D\n,x max ) = 1 + / dxx n - 1 I f (x)D T (n,x). (79) 
Jo 

By this equation the probability for no-branching Pr^(x max ) is that 

Pr^[x max ) = l/D (l,x max ) = exp(- / dxP(x)). (80) 

Jo 

While the probability Pr^{x)dx for the first branching at [x,x + dx](x < x max ) is 
given by 

Pr E (x)dx = dxI f (x)D (l,x)/D (l,x max ) 

= dxl f {x) exp(- / dyl f (y)). (81) 



X 
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By replacing x max by x, we can repeat the use of Eqs.flSOD and (|81|). Then we obtain 
the following iterative algorithm for the shower. 

• step(l) Set x m = 1 and x^ = 1. 

• step(2) Calculate a probability of stop, Pr N (x m ) = exp(— J Xm dyV (y)). 

• step(3) Generate a uniform random number £(0 < £ < 1). If £ is less than the 
probability, go to step (5). 

• step (4) Calculate x that satisfies 

£ = exp(- / dyl f (y)). 
Jo 

Replace Xb by x^x. Then set x m = x and go back to step(2). 

• step(5) Finish generating one event. Then calculate for various n and 
accumulate them. 
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Q 2 = IQ^GeV 2 


Q 2 = 10 6 GeV 2 


mom 


Analytic 


Showers 


Analytic 


Showers 


2 


0.96536 


0.96539 ±0.16E-05 


0.95836 


0.95833 ±0.86E-06 


3 


0.94740 


0.94745 ±0.20E-05 


0.93670 


0.93669 ±0.11E-05 


4 


0.93525 


0.93531 ±0.21E-05 


0.92200 


0.92201 ±0.13E-05 


5 


0.92609 


0.92617 ±0.22E-05 


0.91091 


0.91093 ±0.14E-05 


6 


0.91878 


0.91886 ±0.23E-05 


0.90202 


0.90205 ±0.16E-05 


7 


0.91270 


0.91279 ±0.24E-05 


0.89464 


0.89466 ±0.17E-05 


8 


0.90751 


0.90761 ±0.25E-05 


0.88833 


0.88836 ±0.18E-05 


9 


0.90301 


0.90311 ±0.26E-05 


0.88283 


0.88287 ±0.19E-05 


10 


0.89902 


0.89913 ±0.27E-05 


0.87797 


0.87801 ±0.20E-05 


11 


0.89546 


0.89558 ±0.28E-05 


0.87363 


0.87366 ±0.20E-05 


12 


0.89225 


0.89237 ±0.29E-05 


0.86969 


0.86973 ±0.21E-05 


13 


0.88932 


0.88944 ±0.29E-05 


0.86611 


0.86615 ±0.22E-05 


14 


0.88663 


0.88676 ±0.30E-05 


0.86282 


0.86286 ±0.22E-05 


15 


0.88415 


0.88429 ±0.31E-05 


0.85978 


0.85982 ±0.23E-05 


16 


0.88186 


0.88199 ±0.31E-05 


0.85696 


0.85700 ±0.23E-05 


17 


0.87971 


0.87985 ±0.32E-05 


0.85433 


0.85437 ±0.24E-05 


18 


0.87771 


0.87785 ±0.32E-05 


0.85187 


0.85191 ±0.24E-05 


19 


0.87583 


0.87597 ±0.33E-05 


0.84956 


0.84959 ±0.24E-05 


20 


0.87405 


0.87420 ±0.33E-05 


0.84738 


0.84741 ±0.25E-05 


30 


0.86047 


0.86065 ±0.35E-05 


0.83061 


0.83063 ±0.27E-05 


40 


0.85130 


0.85150 ±0.35E-05 


0.81923 


0.81925 ±0.29E-05 


50 


0.84445 


0.84467 ±0.35E-05 


0.81070 


0.81071 ±0.30E-05 


60 


0.83904 


0.83926 ±0.34E-05 


0.80392 


0.80392 ±0.31E-05 


70 


0.83458 


0.83481 ±0.33E-05 


0.79832 


0.79832 ±0.32E-05 


80 


0.83080 


0.83105 ±0.32E-05 


0.79357 


0.79356 ±0.33E-05 


90 


0.82754 


0.82780 ±0.32E-05 


0.78945 


0.78944 ±0.34E-05 


100 


0.82468 


0.82494 ±0.31E-05 


0.78583 


0.78582 ±0.34E-05 



Table 1: Numerical results of analytic calculations and the shower for the effective LL order 
P-function. Analytic results are given by Eq.(38). Here [i 2 = 0.25 x 10~ 6 GeV, and ao = 1/137. 
The same values for them are used in coming tables. 
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mom 


#6i(l - h)x b2 (l - t 2 ) 


(xn - h)(x b2 - t 2 ) 


q 2 /S 


2 


0.93212 ±0.34E-04 


0.93118 ±0.34E-04 


0.93112 ±0.34E-04 


3 


0.89787 ±0.48E-04 


0.89725 ±0.48E-04 


0.89721 ±0.48E-04 


4 


0.87504 ±0.58E-04 


0.87458 ±0.58E-04 


0.87455 ±0.58E-04 


5 


0.85804 ±0.66E-04 


0.85767 ±0.66E-04 


0.85765 ±0.66E-04 


6 


0.84457 ±0.71E-04 


0.84426 ±0.71E-04 


0.84424 ±0.71E-04 


7 


0.83346 ±0.75E-04 


0.83320 ±0.75E-04 


0.83318 ±0.75E-04 


8 


0.82405 ±0.79E-04 


0.82382 ±0.79E-04 


0.82381 ±0.79E-04 


9 


0.81590 ±0.81E-04 


0.81570 ±0.81E-04 


0.81569 ±0.81E-04 


10 


0.80874 ±0.84E-04 


0.80856 ±0.84E-04 


0.80855 ±0.84E-04 


11 


0.80236 ±0.85E-04 


0.80220 ±0.86E-04 


0.80219 ±0.86E-04 


12 


0.79663 ±0.87E-04 


0.79648 ±0.87E-04 


0.79647 ±0.87E-04 


13 


0.79142 ±0.88E-04 


0.79128 ±0.88E-04 


0.79127 ±0.88E-04 


14 


0.78666 ±0.90E-04 


0.78653 ±0.90E-04 


0.78652 ±0.90E-04 


15 


0.78228 ±0.91E-04 


0.78216 ±0.91E-04 


0.78215 ±0.91E-04 


16 


0.77822 ±0.92E-04 


0.77811 ±0.92E-04 


0.77811 ±0.92E-04 


17 


0.77446 ±0.92E-04 


0.77435 ±0.92E-04 


0.77435 ±0.92E-04 


18 


0.77094 ±0.93E-04 


0.77084 ±0.93E-04 


0.77084 ±0.93E-04 


19 


0.76765 ±0.94E-04 


0.76755 ±0.94E-04 


0.76755 ±0.94E-04 


20 


0.76455 ±0.95E-04 


0.76446 ±0.95E-04 


0.76446 ±0.95E-04 


30 


0.74103 ±0.99E-04 


0.74097 ±0.99E-04 


0.74096 ±0.99E-04 


40 


0.72535 ±0.10E-03 


0.72531 ±0.10E-03 


0.72531 ±0.10E-03 


50 


0.71376 ±0.11E-03 


0.71373 ±0.11E-03 


0.71373 ±0.11E-03 


60 


0.70466 ±0.11E-03 


0.70463 ±0.11E-03 


0.70462 ±0.11E-03 


70 


0.69720 ±0.11E-03 


0.69718 ±0.11E-03 


0.69718 ±0.11E-03 


80 


0.69092 ±0.11E-03 


0.69090 ±0.11E-03 


0.69090 ±0.11E-03 


90 


0.68553 ±0.11E-03 


0.68551 ±0.11E-03 


0.68551 ±0.11E-03 


100 


0.68080 ±0.11E-03 


0.68079 ±0.11E-03 


0.68078 ±0.11E-03 



Table 2: Numerical results of the shower and the generator in the effective LL order. Columns 
with Xf,i{l — ti)x b2 (l — t 2 ) and {x b i — t\)(xb2 — t 2 ) are results of the shower for these variables, 
while the third columns are results of the generator for q 2 /S. 
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mom 


D d (n, Q 2 )/(a /27r) 


D d {n, Q 2 ) 


1 


1.6882 


0.1961E-02 


2 


0.8009 


0.9305E-03 


3 


0.5477 


0.6363E-03 


4 


0.4223 


0.4906E-03 


5 


0.3455 


0.4014E-03 


6 


0.2931 


0.3405E-03 


7 


0.2548 


0.2960E-03 


8 


0.2255 


0.2620E-03 


9 


0.2024 


0.2351E-03 


10 


0.1836 


0.2132E-03 


11 


0.1680 


0.1952E-03 


12 


0.1549 


0.1799E-03 


13 


0.1437 


0.1669E-03 


14 


0.1340 


0.1556E-03 


15 


0.1255 


0.1458E-03 


16 


0.1181 


0.1372E-03 


17 


0.1115 


0.1295E-03 


18 


0.1056 


0.1226E-03 


19 


0.1003 


0.1165E-03 


20 


0.0955 


0.1109E-03 


30 


0.0646 


0.7501E-04 


40 


0.0488 


0.5670E-04 


50 


0.0392 


0.4556E-04 


60 


0.0328 


0.3809E-04 


70 


0.0282 


0.3272E-04 


80 


0.0247 


0.2869E-04 


90 


0.0220 


0.2553E-04 


100 


0.0198 


0.2300E-04 



Table 3: Numerical results on Dd(n, Q 2 ) of Eq.(57) in Sect. 8. They are shown in the third column, 
while values in the second one are calculated by dividing them by ao/27r. While in other tables 
the moment starts from 2, results for moment 1 are included here. 



30 



mom 


D'(n)- 1 


Df> A {n) 


{jl) shower 


2 


0.99714 


0.99701 


0.99700±0.82E-05 


3 


0.99450 


0.99443 


0.99442±0.15E-04 


4 


0.99221 


0.99217 


0.99216±0.20E-04 


5 


0.99017 


0.99015 


0.99015±0.24E-04 


6 


0.98835 


0.98834 


0.98833±0.27E-04 


7 


0.98669 


0.98668 


0.98667±0.30E-04 


8 


0.98516 


0.98516 


0.98515±0.32E-04 


9 


0.98375 


0.98375 


0.98374±0.35E-04 


10 


0.98243 


0.98244 


0.98243±0.37E-04 


11 


0.98120 


0.98120 


0.98119±0.38E-04 


12 


0.98004 


0.98004 


0.98003±0.40E-04 


13 


0.97894 


0.97894 


0.97893±0.42E-04 


14 


0.97790 


0.97790 


0.97789±0.43E-04 


15 


0.97690 


0.97691 


0.97690±0.44E-04 


16 


0.97596 


0.97596 


0.97595±0.45E-04 


17 


0.97505 


0.97506 


0.97505±0.46E-04 


18 


0.97418 


0.97419 


0.97418±0.48E-04 


19 


0.97334 


0.97335 


0.97334±0.49E-04 


20 


0.97254 


0.97254 


0.97253±0.49E-04 


30 


0.96577 


0.96578 


0.96577±0.57E-04 


40 


0.96054 


0.96055 


0.96054±0.62E-04 


50 


0.95624 


0.95626 


0.95625±0.66E-04 


60 


0.95258 


0.95260 


0.95258±0.68E-04 


70 


0.94938 


0.94941 


0.94939±0.71E-04 


80 


0.94653 


0.94656 


0.94654±0.73E-04 


90 


0.94395 


0.94399 


0.94397±0.75E-04 


100 


0.94161 


0.94165 


0.94162±0.76E-04 



Table 4: Numerical results on the compensation for the Q 2 -independcnt contributions. The 
second figures show (n)^ 1 , while the third ones are approximated contributions by D^ A {n). 
Results by the shower algorithm are given in the last column 
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mom 


LL ordcrfanalyt.) 


NLL order (analyt.) 


xm (1 — ti )xh2(l — to) 


q 2 /S 


2 


0.92649 


0.92677 


0.92655±0.37E-04 


0.92554±0.38E-04 


3 


0.88754 


0.88799 


0.88789±0.51E-04 


0.88723±0.51E-04 


4 


0.86085 


0.86143 


0.86140±0.57E-04 


0.86090±0.57E-04 


5 


0.84055 


0.84124 


0.84124±0.61E-04 


0.84085±0.61E-04 


6 


0.82421 


0.82498 


0.82500±0.63E-04 


0.82467±0.63E-04 


7 


0.81056 


0.81141 


0.81143±0.65E-04 


0.81115±0.65E-04 


8 


0.79886 


0.79977 


0.79979±0.67E-04 


0.79955±0.67E-04 


9 


0.78863 


0.78960 


0.78963±0.68E-04 


0.78941±0.68E-04 


10 


0.77956 


0.78058 


0.78060±0.69E-04 


0.78041±0.69E-04 


11 


0.77142 


0.77249 


0.77251±0.70E-04 


0.77233±0.70E-04 


12 


0.76404 


0.76516 


0.76517±0.71E-04 


0.76501±0.71E-04 


13 


0.75730 


0.75846 


0.75847±0.72E-04 


0.75832±0.72E-04 


14 


0.75110 


0.75229 


0.75230±0.73E-04 


0.75217±0.73E-04 


15 


0.74536 


0.74659 


0.74660±0.74E-04 


0.74647±0.74E-04 


16 


0.74002 


0.74128 


0.74129±0.74E-04 


0.74117±0.75E-04 


17 


0.73504 


0.73633 


0.73633±0.75E-04 


0.73622±0.75E-04 


18 


0.73036 


0.73168 


0.73168±0.76E-04 


0.73158±0.76E-04 


19 


0.72596 


0.72731 


0.72730±0.77E-04 


0.72721±0.77E-04 


20 


0.72180 


0.72318 


0.72317±0.77E-04 


0.72308±0.77E-04 


30 


0.68965 


0.69123 


0.69120±0.81E-04 


0.69114±0.81E-04 


40 


0.66759 


0.66932 


0.66928±0.83E-04 


0.66924±0.83E-04 


50 


0.65093 


0.65277 


0.65271±0.83E-04 


0.65268±0.83E-04 


60 


0.63759 


0.63952 


0.63946±0.84E-04 


0.63943±0.84E-04 


70 


0.62652 


0.62853 


0.62846±0.84E-04 


0.62843±0.84E-04 


80 


0.61707 


0.61914 


0.61907±0.85E-04 


0.61905±0.85E-04 


90 


0.60885 


0.61098 


0.61090±0.85E-04 


0.61089±0.85E-04 


100 


0.60158 


0.60377 


0.60369±0.86E-04 


0.60367±0.86E-04 



Tabic 5: Numerical results in the NLL order. The second figures show moments in the LL 
order, while the third ones are the NLL order results. Both are calculated analytically. Results 
on Xbi(l — £i)xb2(l — ^2) and q 2 /S obtained by our generator are given in the last two columns. 
S = 10 4 GeV 2 , = 0.25 x 10- 6 GeV 2 . 
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mom 


LL order(analyt.) 


NLL order(analyt.) 


Xbi(l - h)x b2 (l - t 2 ) 


q 2 /S 


2 


0.91309 


0.91341 


0.91318 ±0.53E-04 


0.91216 ±0.55E-04 


3 


0.86757 


0.86809 


0.86795 ±0.66E-04 


0.86729 ±0.67E-04 


4 


0.83658 


0.83726 


0.83717 ±0.74E-04 


0.83668 ±0.74E-04 


5 


0.81315 


0.81394 


0.81387 ±0.79E-04 


0.81348 ±0.80E-04 


6 


0.79436 


0.79525 


0.79518 ±0.84E-04 


0.79486 ±0.84E-04 


7 


0.77872 


0.77969 


0.77962 ±0.87E-04 


0.77934 ±0.87E-04 


8 


0.76535 


0.76639 


0.76632 ±0.90E-04 


0.76608 ±0.90E-04 


9 


0.75370 


0.75481 


0.75472 ±0.93E-04 


0.75451 ±0.93E-04 


10 


0.74338 


0.74455 


0.74446 ±0.95E-04 


0.74427 ±0.95E-04 


11 


0.73415 


0.73537 


0.73527 ±0.97E-04 


0.73510 ±0.97E-04 


12 


0.72579 


0.72706 


0.72696 ±0.98E-04 


0.72681 ±0.98E-04 


13 


0.71818 


0.71949 


0.71938 ±0.10E-03 


0.71924 ±0.10E-03 


14 


0.71118 


0.71253 


0.71242 ±0.10E-03 


0.71229 ±0.10E-03 


15 


0.70471 


0.70610 


0.70598 ±0.10E-03 


0.70586 ±0.10E-03 


16 


0.69871 


0.70013 


0.70001 ±0.10E-03 


0.69990 ±0.10E-03 


17 


0.69311 


0.69456 


0.69444 ±0.10E-03 


0.69433 ±0.10E-03 


18 


0.68786 


0.68935 


0.68922 ±0.11E-03 


0.68912 ±0.11E-03 


19 


0.68293 


0.68444 


0.68431 ±0.11E-03 


0.68422 ±0.11E-03 


20 


0.67827 


0.67982 


0.67968 ±0.11E-03 


0.67959 ±0.11E-03 


30 


0.64245 


0.64421 


0.64405 ±0.11E-03 


0.64399 ±0.11E-03 


40 


0.61806 


0.61997 


0.61979 ±0.11E-03 


0.61975 ±0.11E-03 


50 


0.59973 


0.60175 


0.60156 ±0.11E-03 


0.60153 ±0.11E-03 


60 


0.58513 


0.58724 


0.58705 ±0.11E-03 


0.58702 ±0.11E-03 


70 


0.57304 


0.57524 


0.57503 ±0.11E-03 


0.57501 ±0.11E-03 


80 


0.56277 


0.56503 


0.56482 ±0.11E-03 


0.56480 ±0.11E-03 


90 


0.55385 


0.55616 


0.55595 ±0.11E-03 


0.55594 ±0.11E-03 


100 


0.54599 


0.54835 


0.54814 ±0.11E-03 


0.54812 ±0.11E-03 



Table 6: The same as Table 5 except S. Here S = 10 6 GeV 2 . 
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S = 10 4 GeV 2 


S = 10 6 GeV 2 


S = 10 8 GeV 2 


mom 


generator 


shower 


generator 


shower 


generator 


shower 


2 


-0.13218% 


-0.02363% 


-0.13773% 


-0.02573% 


-0.13908% 


-0.02333% 


3 


-0.08570% 


-0.01171% 


-0.09262% 


-0.01624% 


-0.09109% 


-0.01249% 


4 


-0.06095% 


-0.00395% 


-0.06951% 


-0.01075% 


-0.06612% 


-0.00565% 


5 


-0.04660% 


-0.00012% 


-0.05688% 


-0.00872% 


-0.05194% 


-0.00254% 


6 


-0.03770% 


0.00182% 


-0.04917% 


-0.00843% 


-0.04318% 


-0.00130% 


7 


-0.03167% 


0.00271% 


-0.04438% 


-0.00898% 


-0.03751% 


-0.00107% 


8 


-0.02738% 


0.00300% 


-0.04110% 


-0.00979% 


-0.03350% 


-0.00136% 


9 


-0.02432% 


0.00291% 


-0.03882% 


-0.01073% 


-0.03063% 


-0.00180% 


10 


-0.02203% 


0.00269% 


-0.03720% 


-0.01182% 


-0.02859% 


-0.00239% 


11 


-0.02032% 


0.00233% 


-0.03617% 


-0.01278% 


-0.02700% 


-0.00300% 


12 


-0.01895% 


0.00196% 


-0.03521% 


-0.01375% 


-0.02562% 


-0.00347% 


13 


-0.01767% 


0.00158% 


-0.03461% 


-0.01473% 


-0.02447% 


-0.00396% 


14 


-0.01675% 


0.00133% 


-0.03410% 


-0.01558% 


-0.02386% 


-0.00474% 


15 


-0.01607% 


0.00094% 


-0.03385% 


-0.01643% 


-0.02307% 


-0.00524% 


16 


-0.01538% 


0.00054% 


-0.03342% 


-0.01714% 


-0.02269% 


-0.00590% 


17 


-0.01467% 


0.00027% 


-0.03326% 


-0.01785% 


-0.02229% 


-0.00641% 


18 


-0.01435% 


-0.00014% 


-0.03307% 


-0.01857% 


-0.02187% 


-0.00678% 


19 


-0.01389% 


-0.00041% 


-0.03302% 


-0.01914% 


-0.02159% 


-0.00730% 


20 


-0.01369% 


-0.00083% 


-0.03310% 


-0.01986% 


-0.02160% 


-0.00798% 


30 


-0.01288% 


-0.00420% 


-0.03415% 


-0.02530% 


-0.02216% 


-0.01300% 


40 


-0.01300% 


-0.00642% 


-0.03549% 


-0.02871% 


-0.02369% 


-0.01672% 


50 


-0.01348% 


-0.00827% 


-0.03673% 


-0.03141% 


-0.02525% 


-0.01966% 


60 


-0.01407% 


-0.00969% 


-0.03814% 


-0.03355% 


-0.02690% 


-0.02226% 


70 


-0.01448% 


-0.01082% 


-0.03911% 


-0.03529% 


-0.02850% 


-0.02451% 


80 


-0.01486% 


-0.01163% 


-0.04018% 


-0.03681% 


-0.02968% 


-0.02619% 


90 


-0.01522% 


-0.01228% 


-0.04100% 


-0.03794% 


-0.03082% 


-0.02766% 


100 


-0.01557% 


-0.01292% 


-0.04194% 


-0.03921% 


-0.03194% 


-0.02892% 



Table 7: Differences between the analytic results and the generator's or the shower's ones. The 
ratios of them are shown. The columns with q 2 /S denote the ratios by the generator. We ne- 
glect the statistical error. Here the shower and the generator compensate for the <5 2 -independcnt 
contributions. 
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mom 


Eq.(76)/( ao /2^) 2 


Eq.(76) 


1 


3.3366 


0.4503E-05 


2 


4.4792 


0.6045E-05 


3 


5.0860 


0.6864E-05 


4 


5.4672 


0.7378E-05 


5 


5.7301 


0.7733E-05 


6 


5.9227 


0.7993E-05 


7 


6.0702 


0.8192E-05 


8 


6.1868 


0.8350E-05 


9 


6.2813 


0.8477E-05 


10 


6.3594 


0.8583E-05 


11 


6.4252 


0.8671E-05 


12 


6.4812 


0.8747E-05 


13 


6.5296 


0.8812E-05 


14 


6.5718 


0.8869E-05 


15 


6.6089 


0.8919E-05 


16 


6.6418 


0.8964E-05 


17 


6.6712 


0.9003E-05 


18 


6.6976 


0.9039E-05 


19 


6.7214 


0.9071E-05 


20 


6.7430 


0.9100E-05 


30 


6.8840 


0.9291E-05 


40 


6.9574 


0.9390E-05 


50 


7.0023 


0.9450E-05 


60 


7.0327 


0.9491E-05 


70 


7.0546 


0.9521E-05 


80 


7.0711 


0.9543E-05 


90 


7.0840 


0.9560E-05 


100 


7.0944 


0.9574E-05 



Table 8: Numerical results on Diff of Eq.(76) in Appendix C. These differences are shown in 
the third column, while values in the second one are calculated by dividing them by (a /27r) 2 - 
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p l 



P 2 



Fig. 1 

Figure 1: The schematics of the annihilation process. The electron(positron) has the momentum 
Pi{P2) initially, and doespi(p2) after the branching process. The momentum of the virtual photon 
is denoted by q. 
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